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120.  abstract 
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Ideas  and  concepts  about  Cracture  and  crack  propagation  in  overconsolidated 
clays  have  been  part  of  gectechnical  engineering  thinking  for  years^  Those  con¬ 
cepts,  however,  were  never  formalized,  primarily  due  to  the  overwh^ming  hold 
plasticity  theory  had  an  still  h||5»  on  the  profession.  To  the  knowledge  of  the 
authors  this  is  the  first  time"Jhe  mechanical  behavior  of  stiff  clays  is  studied 
and  analyzed  through  the  discipline  of  fracture  mechanics.  ' 
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A  test  Is  developed  to  determine  the  critical  stress  Intensity  factor^K^ 

for  the  opening  mode.  The  specimen  with  the  configuration  of  a  hollow  disk  and  a 
notch  cut  In  Its  Internal  surface  Is  fatigued  under  cyclic  internal  pressure, 
then  failed  to  determine  the  pattern  of  crack  propagation.  The  critical  length 
of  the  crack  Is  •thetF'used  to  determine  for  various  degrees  of  over  consolida¬ 
tion*  Vfhlle  overconsolidation  leads  to  an  Increase  In  strength,  -It-'wa*- 
-thoe^t  resulted  In  a  noticeable  decrease  In  the  value  of 

i.iw  ' 

A  second  test  was  developed  to  determine  the  critical  stress  Intensity 
factor  for  the  shearing  mode.  -l-E — ta  iiotewurthy— frha^'Vhls  Is  the  first  time 

a  non-amnlguous  mode  TI  fracture  has  been  obtained  for  any  material.  The 
specimen  with  the  configuration  of  a  thin  long  hollow  cylinder  with  a  notch  cut 
Into  It,  Is  subiected  to  hydrostatic  pressure  and  pure  torsion.  i  Here  too,  the 
knowledge  of  the  critical  crack  length  and  shearing  stress  allows  Me  to  get 

for  a  large  spectrum  of  oyerconsolldatlon  ratios.  Here,  however,  nfhe  Influence 

I  of  overconsolidation  on  Kil*_  was  much  smaller  than  that  observed  on  K^. 

lie  1C  ^ 

Most  Intriguing  are  the  observed  damage  zones  which  have  very  little  1' 
anything  to  do  with  the  predictions  of  plasticity  theory.  A  honey  combed  network 
of  fissures  was  observed  enclosing  undamaged  material  with  localized  discon¬ 
tinuities  emanating  from  the  tip  of  the  crack.  Here  damage  zone  theory  rather 
than  plasticity  Is  the  most  promising  tool  In  the  process  of  understanding  the 
deformation  and  failure  of  fissured  clays. 

j— 1 

^IC  ^IIC  catastrophic  failure  and  do  not  describe  the 

phenomenon  of  slow  crack  propagation.  The  klnetlc.8  of  crack  propagation  for  mode 
*I  was  studied  by  direct  surface  observation^ "  Those  observations  were^compared  tc 
crack  lengths  calculated  from  crack  opening  displacement  measurements.,  The  rela¬ 
tion  between  the  rate  of  crack  propagation  with  respect  to  the  number  of  cycles, 
and  the  energy  release  rate  which  was  deduced.  Is  at  the  core  of  the  development 
of  constitutive  equations  for  crack  propagation;  provided  it  Is  verified  ir 
further  experiments  on  different  clays. 

A  thermodynamic  criterion  of  stability  Is  applied  to  the  stability  of  at 
Infinite  slope  with  a  crack  parallel  to  the  surface. 
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PART  I:  IRTRODOCTION 


Ta  -  Soil  Mechanics  and  Fracture  Mechanics  -  A  Review  of  Some  Relevant  Work 

Cracks  are  very  often  observed  In  stiff  clays.  Some  appear  during  a 
shrinkage  process;  some  appear  during  an  expansion  process;  and  some  are 
Initiated  by  high  levels  of  stress  concentration  at  low  nominal  stress 
levels.  Once  a  crack  has  been  Initiated  subsequent  propagation  may  follow 
any  of  the  three  modes  shown  In  Fig.  1.  Mode  I,  the  opening  mode,  is  of 
Importance  In  hydraulic  fracturing  and  modes  II  and  III,  the  shearing 
modes,  are  of  Interest  In  all  problems  Involving  sliding. 

Ideas  and  concepts  about  fracture  and  crack  propagation  have  been  In 
the  back  of  the  mind  of  geotechnical  engineering  researchers  for  years. 
Often  the  terminology  is  different  from  that  used  by  classical  fracture 
mechanicians,  but  the  problems  are  common.  The  geotechnical  engineering 
profession  Is  quite  familiar  with  stiff  fissured  or  jointed  clays  where  the 
discontinuities.  In  contrast  to  the  Griffith  cracks,  are  found  on  a 
macroscopic  scale  (Fig.  2).  The  possible  geological  processes  that  may 
lead  to  the  formation  of  fissures  has  been  discussed  by  Fookes  [7],  A 
preliminary  classification  of  different  types  of  discontinuities  which 
constitute  planes  of  weakness  has  been  given  by  Skempton  and  Petley  [15] 
and  Skemptlon,  Schuster  and  Petley  [16].  The  importance  of  size  effect  on 
the  undralned  strength  of  stiff  fissured  clay  was  demonstrated  by  Bishop 
and  Little  [2]  who  showed  that  the  strength  of  London  clay  determined  by 
In-sltu  shear  box  testing  on  a  2  ft.  by  2  ft.  square  sample  is  only  55 
percent  of  that  given  by  laboratory  unconsolidated  undrained  trlaxlal  tests 
on  1.5  In.  diameter  samples.  Similar  strength  decrease  was  found  for  blue 


Vertical  Section  Through  Typical  Fissured  Clay 


London  clay  by  Simons  [13].  For  a  glaclo-lacustrlne  clay  north  of  lake 
Erie,  Lo,  adams  and  Seychuk  Il9]  found  that  the  size  effect  was  more 
pronounced.  The  phenomenon  is  therefore,  not  only  a  mathematical  prob¬ 
ability,  but  a  physical  reality  of  no  Infrequent  occurence.  The  remarkable 
work  of  Lo  [9]  on  the  operational  strength  of  fissured  clay  has  in  it  most 
of  the  elements  on  which  classical  fracture  mechanics  rests. 

A  characteristics  common  to  "most"  stiff  clays  is  that  they  are  over¬ 
consolidated  and  therefore  they  have  the  capacity  of  decreasing  their 
shearing  strength  with  increasing  deformation.  Deformation  is  time 
dependent  and,  as  pointed  out  by  Morgenstern  [10],  one  has  to  distinguish 
betwen  two  aspects  of  the  time  dependent  behavior  of  soils.  These  are 
creep  and  progressive  failure.  A  soil  exhibiting  creep  at  a  given  stress 
level  will  still  be  able  to  sustain  a  higher  stress  level  although  possibly 
with  a  higher  strain  rate.  The  best  documented  evidence  for  creep  in  soils 
is  secondary  consolidation  and  field  evidence  shows  that  secondary 
consolidation  in  stiff  clays  can  be  substantial  and  depend  upon  the  history 
of  loading  [A].  Progressive  failure  relates  to  the  fact  that  the  available 
shearing  resistance  of  a  soil  mass  can  be  between  the  peak  and  the  residual 
value.  The  strength  can  be  attacked  in  various  ways  giving  rise  to  dis¬ 
placements  which  vary  with  time.  When  the  shearing  strength  drops  below 
the  shearing  stress,  catastrophic  movement  can  ensue.  Both  creep  and 
progressive  failure  are  dependent  upon  the  loading  path. 

In  his  Terzaghl  lecture  B.ierrum  [3]  states  that  in  order  to  start  a 
progressive  failure  there  must  be  a  discontinuity  somewhere  in  the  clay 
mass  where  failure  can  be  Initiated  and  the  deformations  required  for 
further  development  can  be  produced.  Stress  concentrations.  In  other  words 
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local  high  shearing  stresses  exceeding  the  peak  are  given  as  a  prerequisite 
for  the  development  of  a  continuous  failure  surface  by  progressive  failure. 
The  more  overconsolidated  the  clay,  the  greater  Its  content  of  recoverable 
strain  energy  and  the  greater  the  danger  of  progressive  failure.  Bishop 
[1]  defined  brittle  clays  by  means  of  a  brittleness  parameter  defined  by 
(Fig.  3) 


which  for  overconsolidated  soils  varies  with  stress.  He  also  suggested 

that  brittleness  may  be  considered  In  terms  of  the  additional  energy 

required  to  make  the  shear  failure  progress  from  the  peak  to  the  residual 

state  (Fig.  3).  This  is  represented  by  the  shaded  area  E^.  If  the  sample 

had  shown  no  drop  In  strength  after  the  peak,  the  work  done  would  have  been 

represented  by  E^.  The  difference  between  these  areas  expressed  as  a  ratio 

of  E  may  he  a  useful  measure  of  brittleness  also.  It  was  called  the 
P 

rupture  index  (Fig.  3); 


E  -  E- 
P  f 


Rupture  Index 


Bishop  cautioned  against  evaluating  any  parameter  based  on  the  shape  of  the 


stress  strain  curve,  partlcularllly  after  the  peak,  because  of  the  tendency 
for  failure  to  develop  In  a  single  thin  shear  zone  In  many  samples  which 
show  marked  dllatancy  or  marked  structural  breakdown  or  shear. 

The  mechanism  of  propagation  of  a  shear  surface  In  a  fissured  clay  may 
well  prove  to  Involve  concepts  similar  to  those  put  forward  by  Griffith. 
These  Include  energy  absorbed  In  forming  the  surface  and  released  by 
relieving  the  stresses  In  the  locally  highly  stressed  material  adjacent  to 
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FIG.  15  Crack  Length  Versus  Nuaber  of  Cycles  for  Four  Tests 
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lie  -  Kinetics  of  Crack  Propagation 

To  study  the  kinetics  of  crack  propagation  four  tests  were  conducted 
during  which  the  crack  length  was  measured  as  the  number  of  internal 
pressure  cycles  Increased.  Fig.  15  shows  the  results  of  those  tests. 
Average  values  are  plotted  in  Fig.  16  together  with  a  measure  of  the  amount 
of  scatter.  As  previously  Indicated  the  crack  length  is  measured  on  the 
face  of  the  disk  with  a  follower  attached  l''  an  LVDT.  There  Is  no 
guarantee  that  this  crack  length  extends  all  the  way  across  the  material. 
Thus,  an  alternate  method  for  recording  the  crack  length  was  sought.  This 
method  Is  based  on  a  simple  relation  between  the  crack  length  and  the  crack 
opening  6  at  a  certain  point.  According  to  Castlgllano’ s  second  theorem, 
for  a  linearly  elastic  material,  the  displacement  6  caused  by  a  load  Q  Is 
given  by 


W  _  ^ 

3Q  "  “  3Q 


(6) 


where  U  Is  the  strain  energy  of  the  solid  and 
Let  us  consider  the  disk  shown  In  Fig 
can  be  expressed  as: 

n  -  n 


no  crack 


9n 

3Q 


n  the  potential  energy. 

17.  The  potential  energy 

(7) 


The  derivative  -  —  is  the  energy  release  rate  J,  for  linear  elasticity  and 

(7  Jl  X 

can  be  expressed  in  terms  of  the  stress  Intensity  factor: 


3il  E  *^1 


(8) 


To  compute  the  crack  opening  6  at  a  distance  c  from  the  Inner  face  of  the 
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(10^ 


Green's  function  G^(x)  for  the  hollow  disk  with  a  crack  emanating  from  the 
hole  la  not  known.  As  a  first  approximation  use  was  made  of  Green's 
function  for  the  stress  Intensity  factor  In  the  strip  shown  In  Fig.  12: 


“k  -  ii  '<!•  -t’ 

X  iL 

where  F(-,  -)  Is  a  geometric  correction  factor  given  In  reference  17  page 
I  w 

l~n .  To  apply  Eq.  3  to  the  hollow  disk  we  set  (b-a)  equal  to  M.  b  and  a 
are  the  outer  and  inner  radii  of  the  disk.  The  elastic  stress  distribution 
along  a  radius  In  a  hollow  disk  subjected  to  Internal  pressure  is  given  by 


(Fig. 13). 


oa2  Vi  2 

OgCa  +  x)  »  — 2 - 2  {1  +  ( - )  ] 

b-a  a  +  X 


Using  the  superposition  principle  and  equations  2,  3  and  4,  the  stress 
Intensity  for  a  hollow  disk  with  a  crack  emanating  from  the  hole  Is 


K  = 

^  /ttI  (b^  -  a?- 


X 


Numerical  Integration  of  Eq.  5  leads  to  the  curve  shown  in  Fig.  12. 
If  a  and  p  are  replaced  by  their  critical  value  and  and  substituted 
In  Eq.  5,  the  critical  value  K  of  the  stress  Intensity  factor  Is 

.L  \f 

obtained.  For  the  tests  conducted  In  this  investigation,  this  value  as  a 
function  of  the  overconsolidation  ratio  Is  shown  In  Fig.  lA.  Is  a 

basic  parameter  to  be  used  In  design  when  failures  such  as  those  occurring 
during  hydraulic  fracturing  are  studied. 

Damage  zone  propagation  studies  are  presently  In  progress  to  clarify 
the  dependence  of  K  on  the  preconsolidation  pressure. 

^  V« 
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Kt/K 


FIG.  11  Calculation  of 


J/(b-a) 

FIG.  12  Datanlnation  of  for  a  Hollow  Dlak 
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overconsolidation  ratio  were  mounted  In  the  testing  fixture  as  shown  In 
Fig.  8  and  sinusoidally  varying  air  pressure  sent  through  the  air  Inlet. 
Under  cyclic  pressure  the  crack  propagates  In  a  stable  manner  as  long  as 
the  pressure  Is  moderate  and  the  crack  length  below  the  critical .  The 
crack  was  allowed  to  propagate  half  way  through  thickness  to  allow  for  the 
study  of  damage  zones.  Crack  propagation  was  followed  with  a  pointer 
attached  to  an  LVDT.  The  crack  opening  displacement  (COD)  was  measured 
with  a  transducer  mounted  on  the  specimen.  The  cyclic  loading  was  then 
stopped  and  a  sudden  high  pressure  was  applied  Inside  the  membrane  which 
caused  failure  In  a  very  short  time.  An  X-Y  plotter  gave  a  graph  of  the 
Inside  pressure  versus  the  COD.  The  maximum  pressure  Is  the  critical 
pressure.  The  critical  length  was  obtained  from  the  examination  of  the 
fracture  surface. 

Fig.  9  shows  the  whole  apparatus  designed  and  built  for  this  purpose. 
Fig.  10  shows  the  details  of  the  micrometer  driven  pointer  following  the 
crack  propagation. 

Ilb  -  Calculation  of 

The  most  general  and  elegant  way  to  calculate  is  through  the  use  of 
Green's  function.  Let  Green's  function  for  the  stress  Intensity 

factor  for  a  given  sample  geometry;  l.e.,  the  expression  for  the  stress 
Intensity  factor  caused  by  the  unit  double  forces  applied  at  the  point  x 
(Fig.  11); 

Kf  -  Gj^(x)  (1) 

Then  for  an  arbitrarily  distributed  traction  o(x)  normal  to  the  crack's 

edges,  the  stress  Intensity  factor  Is  given  by, 

l+z 

Ki  ■  J  Gj^(x)  o(x)dx 


(2) 


FIG.  10  Crack  Propagation  mad  COD  NeaaureBents 


FIG.  7  Specimen  for  Mode  I  Test 


FIG.  8  Mode  I  Test  Fixture 


PART  II:  DETERMIIIATION  OF  THE  CRITICAL  TALDB  OP  THE  STRESS  IllTERSITT 


FACTOR  FOR  MODE  I  FRACTURE 


Fracture  Modes  I  and  II  are  respectively  distinguished  by  symmetry  and 
antisymmetry  of  stress,  strain  and  displacement  fields  with  respect  to  the 
crack  plane.  For  Mode  I  fracture  the  sample  configuration  used  Is  that  of  a 
hollow  disk  with  a  precut  notch  or  stress  concentration  cut  In  the  Inside 
edge.  Static  or  dynamic  pressure  Is  applied  to  the  Inner  surface  through  a 
rubber  membrane. 

For  linear  elastic  media,  the  stress  Intensity  factor  is  uniquely 
defined  by  the  stress  and  the  length.  Thus,  Its  critical  value  can  be 
obtained  once  the  value  of  the  critical  stress  and  that  of  the  critical 
crack  length  at  which  unstable  fracture  takes  place  are  known. 

Ila  -  Sample  Preparation  and  Testing  Equipment 

Kaollnlte  clay  with  liquid  and  plastic  limits  of  36.3  percent  and  37.5 
percent  respectively  was  mixed  with  dealred  distilled  water  at  twice  the 
liquid  limit.  The  slurry  was  consolidated  in  a  large  metal  consolldometer 
under  a  series  of  pressures  namely  414,  517,  620  and  723  kPa  to  provide  a 
range  of  overconsolidation  ratios.  After  equilibrium  was  reached  the  clay 
was  allowed  to  rebound  under  atmospheric  pressure  leading  to  overconsollda- 
tlon  ratios  of  4,  5.1,  6.1  and  7,1  respectively.  Fig.  6  shows  the  consoll¬ 
dometer  set-up. 

The  blocks  of  clay  were  sliced  horizontally  to  yield  specimens  25mm 
thick.  A  simple  rig  was  used  to  punch  a  35.6mra  hole  at  the  center  of  the 
discs  and  to  cut  a  radial  notch  10mm  deep  (Fig.  7).  They  were  covered  with 
silicone  oil  wrapped  and  stored  In  a  fog  room.  The  specimens  with  a  given 
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their  density  and  their  statistical  distribution  In  the  material 
Ic  -  Goals  of  the  Present  Study 


Following  a  path  slmlllar  to  that  taken  for  other  materials,  this 
study  seeks  to  develop  standard  tests  for  the  determination  of  and 

for  stiff  clays.  Zones  of  microcracking  are  studied  as  well  as  the 
kinetics  of  crack  propagation.  This  and  future  Investigations  well 
hopefully  lead  to  realistic  constitutive  relations  and  lay  a  foundation  for 
better  understanding  the  behavior  of  stiff  clays  under  static  and  dynamic 
loading. 


very  small  peak.  High  stress  concentrations  or  large  strains  are  needed  to 
Initiate  cracks  In  such  clays.  Often,  however,  inclusions  such  as  coarse 
particles  from  stress  concentrators  for  cracks  to  start  and  propagate. 

Due  to  the  presence  of  a  crack  the  strength  of  the  soil  Is  decreased. 
As  the  size  of  the  crack  Increases  with  time  additional  degradation  of  the 
strength  takes  place  (Fig.  5).  After  a  certain  time  It  becomes  so  low  that 
the  soil  cannot  withstand  accidental  high  loads  that  may  occur  In  service. 
From  this  moment  on  the  soil  Is  liable  to  fail.  On  the  other  hand,  If 
accidental  high  loads  do  not  occur,  cracks  may  continue  to  grow  with  time 
until  the  strength  falls  below  the  service  load.  With  respect  to  Fig.  5, 
fracture  mechanics  should  help  answer  the  following  questions: 

1  -  What  Is  the  strength  as  a  function  of  the  crack  size? 

2  -  What  size  of  crack  can  be  tolerated  at  the  expected 

service  stresses  (static  or  dyanmlc)?  And  what  Is  the 
critical  crack  size? 

3  -  Under  a  given  state  of  loading,  how  long  does  it  take 

for  a  crack  to  grow  from  a  certain  Initial  size  to  the 
critical  size? 

4  -  How  do  cracks  coalesce  and  ultimately  reduce  the  strength 

to  its  residual  value? 

Many  theories  which  have  been  found  satisfactory  for  both  brittle  and 
ductile  materials  have  never  been  tried  for  clays.  The  aim  of  this  research 
is  to  examine  the  mechanical  behavior  of  clay  from  the  point  of  view  of 
fracture  mechanics.  Since  failure  generally  appears  In  the  form  of  a 
surface,  one  has  to  start  by  studying  the  single  crack  and  Its  propagation, 
followed  by  the  coalescence  of  cracks;  a  topic  which  obviously  Involves 
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presently  the  subject  of  some  research  and  may  become  another  useful  tool 
In  the  geotechnical  engineering  arsenal.  Some  results  suggest  that  the 
mechanism  Involved  Is  that  of  a  tensile  failure  resulting  from  water 
pressure  acting  on  a  crack  In  a  soil  medium.  A  discontinuity  Is  needed  and 
the  soil  at  the  tip  of  the  crack  deforms  plastically. 

More  recently  researchers  have  turned  their  attention  to  the  behavior 
of  clays  under  large  cyclic  loadings.  Under  such  loadings  cracks  can  be 
Initiated  as  a  result  of  plastic  deformations.  Even  if  the  nominal 
stresses  are  well  below  the  yield  limit,  locally  the  stresses  may  be  above 
yield  due  to  stress  concentrations  such  as  those  occurring  In  simple  shear 
devices.  Cracks  can  coalesce  after  a  number  of  cycles  leading  to  fatigue 
fracture. 

lb  -  Some  Questions  to  be  Answered  by  Fracture  Mechanics 

Under  stress  clays  seem  to  have  in  them  the  capability  of  responding 
as  elastic,  viscous  or  plastic  materials.  It  all  depends  on  the  level  and 
the  rate  at  which  the  stresses  are  applied.  They  can  be  In  a  brittle  state 
or  In  a  plastic  state;  they  can  be  Intact  or  fissured;  they  can  be 
Insensitive;  they  can  be  strain  hardening  or  strain  softening.  High 
overconsolldatlon  ratios  lead  to  clays  which  develop  fissures  upon  rebound, 
and  such  fissures  can  be  visible  to  the  naktd  eye;  they  have  a  residual 
strength  much  lower  than  the  peak.  Moderate  overconsolldatlon  ratios 
result  In  clays  which  "may"  have  In  them  flaws  or  micro-cracks  Invisible  to 
the  naked  eye,  but  whose  presence  expresses  itself  through  the  existence  of 
a  peak  and  a  residual  strength;  and  through  an  Influence  of  the  size  of  the 
test  specimen  on  the  strength.  Low  overconsolldatlon  ratios  result  In 
clays  whose  behavior  Is  of  the  strain  hardening  type  or  at  best  presents  a 
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Skempton  [14]  recognized  five  stages  In  the  shearing  deformation  of  a 
stiff  clay.  The  first  one  before  the  peak  Is  one  of  continuous  non- 
homogeneous  strain.  In  the  second  which  occurs  at  or  just  before  the  peak 
Riedel  shears  are  formed  (Fig.  4).  They  lie  en  echelon  at  an  Inclination 
usually  between  10  to  30  degrees  to  the  direction  of  general  movement  (the 
D  axis)  and  conjugate  shears  R*  are  sometimes  seen.  With  further  movement 
a  third  stage  Is  soon  reached  at  which  slip  along  the  Riedel  shear  Is  no 
longer  kinematically  possible  and  the  clay  Is  compelled  to  develop  new  slip 
surfaces  parallel  or  subparallel  to  the  D  axis.  These  are  the  displacement 
shears.  With  greater  movement  the  displacement  shears  extend  and  even¬ 
tually  In  the  fourth  stage  some  of  them  link  up  to  form  a  slip  surface. 
This  surface  Is  undulating  since  the  shears  Involved  were  not  originally 
all  In  line.  Thrust  shears  typically  Inclined  at  about  160=  to  the  D  axis 
also  tend  to  develop  In  the  third  and  fourth  stages.  In  the  fifth  stage 
the  slip  surface  undergoes  appreciable  flattening  as  a  result  of  still 
greater  movements.  Once  a  slip  surface  has  formed  subsequent  movements  are 
mostly  concentrated  along  It;  and  large  displacements  can  then  take  place 
without  a  fundamental  change  In  the  pattern  of  shears,  the  development  of 
shear  Is  accompanied  by  particle  orientation.  Morgenstern  and  Tchalenko 
[11]  and  Yong  and  Warkentln  [18]  have  shown  that  the  slip  surface  consists 
of  a  thin  band  In  which  the  particles  are  strongly  oriented  more  or  less  In 
the  direction  of  the  movement.  A  flattened  slip  surface  must  possess  the 
minimum  shear  resistance  and  this  Is  the  residual  strength  of  the  clay. 

An  Interesting  test  that  Is  sometimes  conducted  In  the  field  In 
conjunction  with  permeability  tests  Is  hydraulic  fracturing.  It  has  been 
used  for  years  to  determine  In-sltu  earth  pressures.  Such  a  test  Is 
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FIG.  17  Derivation  of  Crack  Opening  Displacement 


FIG.  18 
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a  =  1*8  cm 
b=5.5  cm 


o(a+x)  =  (  0.063+3. 36/(a+x)2)p 

The  Use  of  Tada's  Strip  to  Derive  an  Approximate  Solution 


disk,  we  apply  the  fictitious  forces  Q  then  differentiate  n  for  0*0. 


Since  n  ,  does  not  depend  on  Q,  then 

no  crack 


0  suq 


Therefore 


%=0  E  Jo  3Q  0 


(9) 


(10) 


where  Is  the  stress  Intensity  factor  due  to  Internal  pressure  without  Q 
and  Kq  Is  the  stress  Intensity  factor  for  Q  alone. 


Thus , 


6  =  (~fe  f  K  •  dH) 

E  J  P  Q  Q-1 


(11) 


Here  again,  using  Green's  function  for  a  strip  given  by  Tada  ll7]  we 
have  (Fig.  18) 


(12) 


K  =  F(-  -) 

Q  /IFF  H’  w^ 


(13) 


2(1  - 


(lA) 


Eq.  14  was  numerically  Integrated  and  the  results  shown  In  Fig.  19.  Crack 
opening  displacements  (COD)  are  measured  at  the  tip  of  the  notch  which  Is 
10mm  from  the  Inside  edge  of  the  disk.  Fig.  20  shows  the  COD  for  four  tests 
versus  the  number  of  cycles,  and  Fig.  21  the  average  of  those  tests  and  the 
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scatter  Involved.  Combining  figures  19,  20  and  21  one  can  deduce  the 
relation  between  the  crack  length  and  the  number  of  cycles.  This  relation 
is  shown  in  Fig.  22,  together  with  the  experimentally  measured  one.  The 
two  curves  are  slgnflcantly  different.  This  may  be  due  to  the  fact  that 
the  COD  versus  the  crack  length  was  obtained  from  a  strip  Instead  of  from  a 
disk.  Apparently  the  error  due  to  such  an  approximation  Increases  with  the 
crack  length.  Also  linear  elasticity  may  have  limited  applicability  to  the 
material.  Both  these  limitations  can  be  eliminated  by  further  theoretical 
developments.  In  addtlion  the  experimental  length  of  the  crack  was 
obtained  from  outside  observation  on  the  face  of  the  disk.  Cross  sections 
%rould  help  insure  that  this  length  indeed  extends  throughout  the  thickness. 

The  kinetics  of  crack  propagation  in  terms  of  crack  propagation  rate 
versus  crack  length  obtained  from  direct  observation  is  shown  in  Fig.  23. 
The  relation  between  ^  and  I  strongly  depends  on  the  specimen's  geometry 
and  thus  does  not  represent  a  material  behavior.  The  same  time,  the  energy 
release  rate  J  accounting  for  specimen  geometry,  crack  length  and  loading 
conditions  is  commonly  used  as  an  Invariant  parameter.  Fig.  24  shows  the 
energy  release  rate  versus  the  crack  length.  ^  versus  J  is  then  plotted  in 
Fig.  25.  This  last  relation,  if  verified  in  further  experiments,  is  at  the 
core  of  the  development  of  constitutive  equations  for  crack  propagation. 
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Crack  Length,  L(mm) 


FIG.  22  Coaparlson  Between  Neeeured  end  Calculated  Crack  Lengths 
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Crack  Length  Bata,  dl-<IN(mm) 
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FIG.  23  Crack  Length  Rate  Veraua  t 
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PART  III:  DBTERMINATICM  OF  THE  CRITICAL  VALUE  OF  THE  STRESS 


INTENSITY  FACTOR  FOR  MODE  II  FRACTURE 

In  this  mode  the  displacement  of  the  crack  surfaces  Is  parallel  to  the 
plane  of  the  crack.  the  sample  configuration  used  Is  that  of  a  thin  long 
hollow  cylinder  with  the  dimensions  shown  In  Fig.  26.  Once  the  length  of 
the  crack  and  the  critical  shearing  stress  are  known  a  solution  due  to 
Erdogan  and  Ratwanl  [6]  allows  one  to  compute 
Ilia  -  Sample  Preparation  and  Testing  Equipment 

Clay  slurries  consolidated  at  approximately  280  kPa  were  first  cut  In 

a  shape  of  a  hollow  cylinder.  A  special  lathe  was  used  to  reduce  them  to 

the  size  shown  In  Fig.  26.  Two  diametrically  opposite  horizontal  notches 

were  cut  In  the  specimens  and  2  layers  of  teflon  tape  Inserted  In  them  to 

prevent  hydrostatic  pressure  from  closing  and  healing  the  crack.  The 

specimens  were  then  placed  In  a  special  cell  where  they  were  one 

dimensionally  consolidated  under  a  given  cell  pressure.  Three  cell 

pressures  were  used  leading  to  vertical  consolidation  pressures  of  1394, 

1174  and  954  kPa;  K  for  this  clay  being  0.47.  The  samples  were  then 

o 

allowed  to  rebound  at  a  cell  pressure  of  207  kPa  leading  to  vertical 
overconsolldatlon  ratios  of  6.7,  5.7  and  4.6. 

A  system  of  pullles  and  actuator  was  used  to  apply  rapidly  a  very  high 
torque  to  the  piston;  while  transducers  Inside  the  cell  measured  the  actual 
torque  transmitted  to  the  sample,  and  the  angle  of  rotation.  These  values 
were  recorded  on  an  X-Y  plotter.  The  maximum  torque  Is  the  critical  value 
needed  to  compute  KTTr-*  ^^8*  shows  a  series  of  failure  surfaces  In  MODE 
II  for  3  of  the  hollow  cylinders  tested.  Fig.  27  shows  the  torsional  set  up 
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FIG.  26  Failure  of  Hollow  Cylinder  in  Mode 


PIG.  27  Mode  II  Torsional  Testing  System 


th  the  hollow  cylinders  of  clay  inside  pressurized  cells, 
lb  -  Calculation  of  Kjjq 

Knowing  the  geometry  of  the  sample  and  the  crack  as  well  as  the 

ritical  torque,  Erdogan's  [6]  solution  yields  the  value  of  K  .  Fig.  2R 

X  J. 

tiows  the  relation  between  and  the  overconsolidation  ratio  for  the 

aollnlte  tested. 

lie  -  Observations  of  Damaged  Zones  in  the  Vicinity  of  Cracks 

One  of  the  most  intriguing  observations  made  during  the  study  Is 
elated  to  the  changes  in  the  material  micro-structure  caused  by  crack 
ropagatlon.  It  is  well  known  that  right  on  the  fracture  (or  slip)  surface 
he  clay  particles  align  themselves  parallel  to  this  surface.  However, 
leyond  this  exceedingly  thin  layer  made  of  the  thickness  of  a  few 
>artlcles,  a  zone  of  damage  is  present.  To  examine  this  zone  of  damage  a 
Large  portion  of  the  hollow  cylinder  surrounding  the  crack  was  first  vacuum 
Impregnated  with  wax.  Then  a  section  of  about  1  cm  in  height  was  cut 
jerpendlcular  to  the  fracture  surface  with  a  diamond  tipped  waferlng  blade 
[Fig.  29).  When  photographed  in  the  direction  indicated  by  the  arrow.  Fig. 
JO  la  obtained.  The  magnification  in  Fig.  30  is  approximately  10,  which  is 
/ery  small  compared  to  what  is  ususally  obtained  in  the  scanning  electron 
nlcroscope.  The  advantage  of  such  a  magnif icatlon  is  that  one  can  see  the 
overall  picture  of  the  damage  and  the  zone  Involved. 

In  Fig.  30,  a  damage  zone  similar  in  shape  to  the  zone  predicted  by 
plasticity  theory  can  be  seen.  However,  contrary  to  plasticity  which 
Involves  a  field  of  slip  or  Luder  lines,  this  zone  consists  of  a  honey¬ 
combed  network  of  fissures  (dark  regions)  enclosing  undamaged  material 
(light  regions).  the  cellular  structure  of  this  region  is  quite  evident 
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PART  IV:  CONCLDSION 


A  framework  for  the  use  of  fracture  mechanics  In  the  study  of  the 
behavior  of  stiff  clays  has  been  established.  Tests  have  been  proposed  and 
then  conducted  for  the  determination  of  K  _  and  kinetics  of 
crack  propagation  In  Mode  I  has  been  examined.  Experimentally  this  Is  the 
first  time  a  pure  Mode  II  has  been  obtained  for  any  material;  and  the 
observations  made  on  the  microstructure  in  the  damage  zone  seem  to  provide 
a  turning  point  In  the  Interpretation  of  failure  phenomena  In  clays. 

Perhaps  the  most  Important  contribution  of  K.  Terzaghl  to  soil  mech¬ 
anics  was  his  bringing  together  In  one  treatise  many  of  the  branches  of  the 
mechanics  of  continuous  media;  and  casting  them  In  a  form  easily  used  by 
engineers.  Fracture  mechanics  Is  a  relatively  young  science  which  has  Its 
place  In  the  interpretation  of  the  behavior  of  soil  masses.  It  Is  hoped 
that  this  research  will  help  develop  new  testing  and  analysis  techniques 
leading  to  more  rational  design  In  fissured  soils. 


QO  ON 


TABLE  I 


Critical  length  of  cracks  for  various  h  and  1 


\1° 

liS. 

15.00 

20.00 

25.00 

30.00 

35.00 

5.0 

33.13 

20.05 

14.11 

11.04 

9.38  ^ 

6.0 

23.01 

13.92 

9.80 

7.67 

6.51 

8.0 

12.94 

7.83 

5.51 

4.31 

3.66 

10.0 

8.28 

5.01 

3.53 

2.76 

2.34 
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Conventionally  the  energy  release  rate  for  plane  strain  and  MODE  II  Is 

2  2 

written  »  (1  -  v  )K^  E,  where  E  and  v  are  Young's  modulus  and 
Poisson's  ratio  respectively.  Thus, 


2(1  -  v^) 


3K 


K 


II 


II 


3Jl 


>  0 


(23) 


and  condition  (21b)  Is  always  satisfied.  Condition  (22)  can  be  rewritten 
In  the  following  equivalent  form 


*^IIS  ”  ‘^IlC 


(2A) 


where  Is  the  fracture  toughness  coefficient.  It  can  be  obtained  for 

Instance  from  tests  conducted  on  hollow  cylinders  In  MODE  II  fracture. 

Substituting  (22)  Into  (24)  we  obtain  the  expression  for  the  critical  crack 

length  21  .  Thus, 
c 

s^iic 

2£  -  - - - -y  (25) 

c  7r(pgh  sin  21)^ 

The  values  of  for  a  range  of  overconsolidation  ratios  can  be  obtained 

from  Fig.  28.  Table  1  gives  values  of  the  critical  lengths  of  crack  which 

would  cause  catastrophic  failure  for  various  depths  h  and  slope  incllna- 

3  3  3/2 

tlons  1.  The  value  of  Pg  Is  17.3  kN/m  and  that  of  Is  156  x  10  Nm  . 

The  dimensions  are  In  meters.  As  one  can  notice  the  results  seem  to  be 


quite  reasonable. 


sufficient  condition  of  stability  (19)  reduces  to 


3Jj^ 


a 


T 


<  0 

const  — 
const 


(20) 


Therefore,  at  constant  temperture  T,  the  necessary  and  sufficient 
conditions  of  Instability  are 


J 


1 


eff 

*  Y 


(a) 


>  0 


const 


(b) 


(21) 


IVb  -  Application  to  the  Stability  of  Infinite  Slopes 

Consider  an  Infinite  slope  with  a  crack  of  length  2Z  parallel  to  the 
surface  and  at  a  depth  h.  It  will  be  assumed  that  there  Is  no  shearing 
resistance  along  the  crack  and  that  for  saturated  overconsolldated 
undralned  clay,  compressive  stresses  do  not  change  the  shear  response.  A 
MODE  II  fracture  Is  under  consideration.  The  shearing  stress  on  any  plane 
parallel  to  the  surface  at  a  depth  h  Is  given  by  x  »  pgh  sln(l)  cos(l), 
where  1  Is  the  Inclination  of  both  slope  and  crack.  Let  us  now  consider  a 
MODE  II  crack  In  an  Infinite  plane  under  the  uniform  shearing  stress  and 
use  the  principle  of  superposition.  If  K__  Is  the  fracture  toughness  of 

ll09 

the  plane  and  K-to  that  of  the  problem  at  hand  then  K  -  “  because 

XXS  XXb  X  Xoo 

of  the  symmetries  Involved  In  an  Infinite  slope.  Thus, 

■  Kjjg  -  xr/ini  -  ipgh/Tl  8ln(2  1) 


(22) 


appears  as  a  generalized  force  reciprocal  to  the  generalized  flux  i. 

When  the  crack  propagates  spontaneously.  When  (y*Rj^  ~ 

•  •  * 

Jj^)>0,  the  crack  can  prorpagate  only  If  D  iiy  Rj^  -  J^),  l.e..  If  the 

•  A 

dissipation  compensates  for  the  negative  entropy  production  Rj^). 

This  Implies  that  the  crack  propagation  is  controlled  by  dissipative  pro¬ 
cesses  like  Irreversible  deformation,  diffusion,  chemical  reaction,  etc. 
Apparently,  the  critical  condition  for  Instability  (stability)  Is 

Jj  -  y\  -  0  (17) 

which  Is  Griffith's  criterion  If  Y  Rj^  is  taken  as  the  surface  energy.  This 

eff 

quantity  In  the  fracture  mechanics  literature  is  referred  to  as  Y  or  the 
critical  energy  release  rate  While  Eq.  17  is  a  necessary  condition 

for  instability,  the  sufficient  condition  results  from  the  universal 
criterion  of  evolution.  For  a  stable  process,  the  following  inequality 
holds  [12]: 


dt 


dX, 

E  Ju  —  < » 

k  k  dt  - 


where  Is  assumed  to  remain  constant. 
Eqs.  16  and  18  yield 


(18) 


For  an  elastic  crack  D  =  0  and 


dt 


3(Ji  -  YR,) 


^  <  0 

o  •  const  — 


const 


(19) 


where  o  Is  the  stress  tensor.  If  the  damaged  zone  does  not  change  as  the 
crack  length  Increases  one  can  assume  that  Y*^!  “  is  constant.  The 
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PAST  IV:  THBRHODTIIAIIIC  CRITERIA  OF  STABILITY  AND  APPLICATION 


TO  THE  SLOPE  STABILITY  PROBLEM 


IVa  -  Theraodynamlc  Criteria  of  Stability 

For  an  Irreversible  process  the  second  law  of  thermodynamics  states 
that  the  rate  of  entropy  production  is  such  that 

®1  *  k  \  ^^5) 


where  is  a  generalized  flux  and  is  a  generalized  force.  The  rate  of 
entropy  production  for  a  propagating  crack  has  be  written  by  Chudnovsky  (5) 


TSj  -  D  +  -  y*R^)  >  0 


(16) 


where,  T  is  the  absolute  temperature,  D  is  the  portion  of  the  irreversible 
work  spent  on  the  development  of  micro  damage  in  the  vicinity  of  the  crack 

tip,  I  is  the  rate  of  crack  growth;  J,  “  “  ^  is  the  energy  release 

1  uJo 

it 

rate,  H  is  the  potential  energy  also  called  Gibbs  potential;  Y  is  the 

specific  enthalpy  of  damage.  It  can  be  looked  upon  as  the  surface  energy  if 

the  damage  consists  of  microcracks  only.  is  a  value  Integral 

* 

characterizing  the  damaged  zone  around  the  crack  tip  (5].  Vhereas  Y 

represents  the  amount  of  energy  needed  for  the  crack  and  its  surrounding 
damaged  zone  to  advance,  characterizes  the  amount  of  accumulated 
potential  energy  available  for  crack  propagation.  Thus  (y*Rj^  - 
represents  the  energy  barrier  which  must  be  overcome  for  the  crack  to 

advance.  Within  the  framework  of  irreversible  thermodynamics,  -  y*Ri) 
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In  the  enlargements  of  Fig.  31. 

In  addition  to  this  damage  zone  Fig.  30  shows  several  localized  bands 
of  discontinuity  of  slippage.  Enlargements  of  one  of  these  localized 
discontinuities  is  shown  In  Fig.  32.  The  localized  discontinuity  has  a 
cellular  structure  similar  to  the  honeycomb  damage  zone,  but  on  a  much 
smaller  scale.  These  localized  discontinuities  are  Inclined  cracks  which 
have  propagated  from  the  notch  tip  as  shown  In  Fig.  33.  The  horizontal 
crack,  rather  than  the  Inclined  cracks,  led  to  catastrophic  failure. 
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FRAGMENTS  OF  HONEYCOMB  DAMAGE 


FIG.  31  Fragments  of  Honeycomb  Damage 


DAMAGE  PATTERN 
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SCHEMATIC  OF  DAMAGE  CROSS  SECTION 


FIG.  28  R 


for  Overconsolidated  Clay 


lie 
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MODE  II  specimen 

29  Sections  Observed  in  Mode  II  Specimen  Fracture 
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